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ABSTRACT 

The  technique  of  speckle  interferometric  measurement  of  asteroids  is  applied  to 
Vesta  and  Pallas,  yielding  diameters  of  330  ±23  and  673  ±33  km,  respectively.  The  improved 
Vesta  speckle  diameter  is  in  excellent  agreement  with  other  measurements  for  Vesta.  However, 
the  derived  Pallas  diameter  is  slightly  larger  than  other  measurements.  This  is  perhaps 
attributable  to  systematic  errors  arising  as  an  object's  angular  diameter  nears  the  seeing  disk 
size.  It  is  also  pointed  out  that  it  is  necessary  to  carefully  normalize  and  center  speckle  frames 
before  applying  the  autocorrelation, cross-correlation  subtraction  method  developed  by  Worden 
et  al.  [Icarus  32.  450  (1977)]. 


I.  INTRODUCTION 

The  use  of  speckle  interferometry  to  obtain  high- 
angular-rcsolution  information  has  proven  valuable  since 
first  proposed  by  Labcyric  (1970).  A  summary  of  the 
uses  and  limitations  of  speckle  technique  has  been  given 
by  Worden  (1977).  The  application  of  speckle  techniques 
to  i he  study  of  asteroids  is  appealing  since  it  provides  an 
alternate  means  of  measuring  their  diameter  to  means 
such  as  radiometric  flux  measurements  and  polarization 
observations.  (See  Morrison  1977  for  a  review.)  In  pre¬ 
vious  work  (Worden  el  al.  1977,  Welter  and  Worden 
1978)  a  method  utilizing  speckle  interferometry  to  derive 
diameters  for  low-surfacc-brightncss  objects  such  as 
asteroids  was  developed  and  used  to  derive  an  angular 
diameter  for  Vesta.  It  is  the  purpose  of  this  report  to 
extend  this  method  towards  deriving  a  more  accurate 
diameter  for  Vesta  as  well  as  obtaining  a  diameter  for 
Pallas. 

II.  DATA  AND  MBTUOD 

Speckle  data  consist  of  a  series  of  short-exposure  (f 
5  0.05  s)  photographs  at  large  image  scales.  On  time 
scales  short  enough  to  freeze  turbulent  motions  in  the 
earth's  atmosphere,  a  large  telescope  behaves  as  a  mul¬ 
tiple-aperture  interferometer  with  a  measurable  mod¬ 
ulation  transfer  function  down  to  the  telescope  diffrac¬ 
tion  limit.  This  means  that  angular  details  on  scales 
similar  to  the  finest  theoretical  scale  may  be  observed. 
The  instrument  used  for  these  observations  was  the  Kitt 
Peak  National  Observatory  4-m  telescope  and  spccklq 
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camcra  described  by  Lynds,  Worden,  and  Harvey 
(1976).  The  data  consisted  of  a  set  of  180  35-mm  Tri-X 
exposures  for  Pallas  and  360  exposures  for  Vesta  ob¬ 
tained  on  3  February  1977  UT,  with  exposure  limes  of 
0.03  s.  To  preserve  the  interferometric  character  of  such 
photos  it  is  necessary  to  limit  the  bandpass  of  observa¬ 
tions,  in  this  case  to  AA  =  300  A  al  5500  A.  Although  not 
specifically  required,  point  source  stars  were  also  ob¬ 
served.  A  summary  of  these  observations  is  given  in 
Table  I. 

The  method  used  to  reduce  these  observations  follows 
that  presented  in  Worden  et  al.  (1977).  The  end  product 
of  this  method  is  the  autocorrelation  of  the  diffraction 
limited  image  O(x.y)  given  mathematically  by 

0(Ax.Ay)*0(Ax.Ay) 

»  S  S 0(x.y)0(x  +  Ax,  y  +  Ay)dxdy.  (I) 

From  this  function  angular  diameters  may  be  obtained 
(but  not  images,  since  the  autocorrelation  loses  phase 
information).  As  shown  by  Worden  el  al.  (1977)  the 
desired  autocorrelation  is  derived  by  computing  the  mean 

autocorrelation  of  a  set  of  speckle  images  _/}(/  =  I . 

/V  for  N  speckle  frames), 

(/HA.T.A y)*fi(Ax.Ay)) 

=  /  J  Ax.y)f(x  +  Ax,  v  +  A y)dxdy.  (2) 

_ Table  I.  Observations  and  results. 


Angular 

determination 

Date 

diameter 

Diameter  (Morrison  1977) 

Object 

(UT) 

(arcscc) 

(knt) 

(km) 

Vesta 

3  Feb  0500 

0.470  ±  0.02 

550  ±  23 

538 

i 

Pallas 

3  Feb  0655 

0.73  ±  0.06 

673  ±  55 

608 

48  Gem 

3  Feb  05 15 

Point  source 

i? 

HR  SM 

3  Feb  0705 

Point  source 
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The  mean  cross  correlation  between  this  same  set  of 
speckle  images  is  then  computed,  and  then  subtracted 
from  (2).  As  shown  in  that  paper,  the  result  is  the  desired 
object  autocorrelation  in  ( I ). 

0(.v..i-)*O(.v.>>)  =  (/i(-v.>)  */(-'.>')) 

i(-v.r)>.  (3) 

Several  aspects  of  the  reduction  method  have  been 
learned  which  were  not  apparent  in  the  original  Worden 
tial.  (1977)  paper.  We  note  here  that  it  is  necessary  to 
normalize  each  frame  so  that  the  sum  of  intensities  is  the 
same  from  frame  to  frame.  The  reason  for  this  procedure 
can  be  seen  by  considering  two  speckle  frames  modeled 
by  identical  Gaussians,  one  with  a  total  area  equal  to  I , 
she  second  reduced  by  a  multiplicative  factor  Fol¬ 
lowing  Braccwel!  (1965),  we  observe  that  the  summed 
area  of  the  autocorrelation  is  equal  to  I1  4-  '/i2  =  1-25. 
Similarly,  the  summed  area  of  the  cross  correlations 
between  the  two  frames  would  be  2(1  X  V2)  =  I.  The 
cross  correlation  should  be  identical  to  the  autocorrela¬ 
tion  in  this  ease  for  our  method  to  work.  Consequently, 
it  is  necessary  to  normalize  each  speckle  frame  to  the 
same  value  for  integral  area;  otherwise  the  summed 
autocorrelation  will  always  be  larger  than  the  cross 
correlations,  even  in  the  absence  of  a  real  speckle  signal. 
A  second  problem  arises  from  the  fact  that  each  frame 
is  not  digitized  in  a  perfectly  centered  manner  relative 
to  the  others.  The  autocorrelations  will  always  be  cen¬ 
tered.  However,  if  the  two  frames  comprising  the  cross 
correlation  arc  not  centered  relative  to  each  other  the 
cross  correlation  will  be  offset.  If  each  frame  in  a  given 
set  is  offset  in  a  random  manner  (his  will  mean  that  the 
summed  cross  correlation  from  this  set  will  be  spread  out 
relative  to  the  autocorrelation,  giving  rise  to  a  spurious 
signal  in  the  autocorrelation,  cross-correlation  sub¬ 
traction.  It  is  thus  necessary  to  center  each  speckle  frame 
before  computing  the  correlation  functions.  This  cen¬ 
tering  was  necessary  to  remove  the  negative  signal  at 
large  autocorrelation  distances,  comparable  to  the  size 
of  the  seeing  image  (~0.5*).  We  found  that  centering 
to  about  0.1*  was  adequate  to  eliminate  this  problem. 

In  practice  the  speckle  frames  for  the  asteroids  and 
stars  arc  scanned  using  the  Kilt  Peak  PDS  microdcnsi- 
tomclcr  onto  a  128  X  1 28-clcmcnl  array.  The  scanning 
aperture  was  100  X  I00j<  with  a  step  size  of  1 00  /a.  This 
procedure  results  in  a  pixel  resolution  of  0.02*  assuming 
a  plate  scale  of  (0.198  +  0.008)*/mm  as  measured  by 
II.  McAlister  at  the  Kilt  Peak  National  Observatory. 
This  resolution  undcrsamplcs  the  diffraction  limit 
frequencies  by  almost  a  factor  of  2.  This  was  necessary 
10  minimize  the  size  of  the  arrays  to  be  processed.  For 
frequencies  near  the  diffraction  limit,  this  would  cause 
aliased  spectra  for  objects  near  the  diffraction  limit  in 

•  jilt,  such  as  point  source  objects.  However,  for  large 

•  objects  like  Vesta  and  Pallas,  there  arc  no  frequencies 
4itcar  the  diffraction  limit,  so  aliasing  is  not  a  problem. 
[These  data  arc  then  converted  from  measured  plate 

density  to  relative  intensities  utilizing  a  standard  curve 

iFIS  FA'i  1  • 
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measured  for  Kodak  Tri-X  film,  finch  frame  is  then 
digitally  Fourier  transformed  using  the  Sacramento  Peak 
Xerox  Sigma  5  computer.  As  discussed  by  Braccwell 
(1965),  the  auto  and  cross  correlations  can  be  computed 
from  the  transforms  F,(u.v)  of  the  frames  fi(x.y)  as 
follows  for  the  autocorrelation: 

/,*/.  =  J'J' 

X  F,(u.t')  /«jr(.VM  +  )v)\tJudr.  (4) 

namely  by  inverse  Fourier  transforming  the  power 
spectrum.  The  cross  correlation  may  be  similarly  ob¬ 
tained  from  the  transforms  of  successive  speckle  frames. 
The  correlations  arc  then  summed  for  20-frame  averages 
and  subtracted  as  detailed  in  F.q.  (3).  This  yields  8  esti¬ 
mates  of  the  unperturbed  autocorrelation  for  the  160 
Pallas  frames  and  16  for  the  Vesta  data.  Similarly,  160 
frames  each  for  point  source  objects  48  Gem  and  HR 
3653  were  reduced. 

To  derive  diameter  estimates  from  these  data,  a 
best-fit  estimate  was  computed  for  each  20-frame  av¬ 
erage  between  the  derived  result  and  theoretical  auto¬ 
correlations  computed  for  various  diameter  uniform 
disks.  Due  to  the  possible  aliasing  problem  near  the 
diffraction  limit,  we  did  not  use  the  derived  point  source 
autocorrelation  as  a  comparison  autocorrelation.  How¬ 
ever,  in  oilier  work  we  found  that  the  point  source  au¬ 
tocorrelation  did  not  differ  significantly  from  the  theo¬ 
retical  one.  We  also  found  that  the  diameters  derived  did 
not  differ  for  the  two  methods  provided  we  only  used  the 
higher-intensity  portions  of  the  program  object  auto¬ 
correlations  (relative  intensity  >0.3)  to  derive  diameter 
estimates.  We  have  assumed  no  limb  darkening  for  these 
objects,  which  is  apparently  a  reasonable  assumption 
(Zcllncr  1977).  The  estimates  derived  in  this  manner  and 
the  errors  determined  from  internal  consistency  in  the 
result  arc  given  in  Table  I  along  with  other  measure¬ 
ments  of  asteroid  diameters. 

The  use  of  photographic  data  presents  several  possible 
difficulties.  The  most  serious  of  these  is  the  nonlinearity 
and  nonuniform  sensitivity  of  the  film.  These  effects  may 
result  in  improper  normalization  of  frames,  as  discussed 
previously.  The  result  of  this  error  would  be  to  reduce  the 
cross  correlations  relative  to  the  autocorrelation,  thus 
increasing  the  measured  diameter  of  the  object.  Another 
effect  which  may  be  present  is  scattered  light  in  the 
system.  This  would  also  result  in  a  possible  increase  in 
the  diameter.  A  sample  of  the  reduced  data  is  shown  in 
Fig.  I ,  a  radial  average  of  the  correlation  difference  for 
Pallas,  Vesta,  and  point  source  comparisons.  Both  of  the 
point  source  reductions  showed  sizes  close  to  the  theo¬ 
retical  diffraction-limited  autocorrelation.  However,  the 
low-intensity  portion  of  the  point  source  reduction  docs 
show  sonic  increase  in  size  which  may  perhaps  be  at¬ 
tributable  to  the  effects  outlined  above.  This  would  be 
expected  to  affect  the  Pallas  diameter  more  severely, 
since  it  has  a  size  nearer  to  the  seeing  disk.  Consequently, 
the  speckle  signal  is  lower  and  more  affected  by  sys¬ 
tematic  errors.  Since  the  measured  diameter  is  slightly 
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Fig.  I.  One-dimensional  radial  average  of  autocorrelation  minus 
cross  correlation  for  one  20- fra  me  set  for  Vesta,  Pallas,  and  48  Gem. 
Each  autocorrelation  has  been  normalized  to  unity. 

larger  than  other  measurements,  this  may  in  fact  be  the 
case.  Another  possible  source  of  error  is  due  to  the  fact 
that  the  asteroids  have  a  very  low  surface  brightness,  and 
the  photon  noise  would  effect  each  telescope  resolution 
element  much  more  than  the  much  brighter  (point  by 
point)  star  images.  However,  reduction  of  lOth-magni- 


tude  star  images  from  a  new  digital  camera  we  art! 
testing  indicates  that  this  is  not  a  serious  problem  eve® 
with  as  few  as  80  frames. 

III.  CONCLUSIONS 

We  have  applied  a  correlation  method  for  the  reduc¬ 
tion  of  speckle  interferometry  data  towards  further 
analysis  of  asteroid  speckle  interferometry  data.  Further, 
we  have  derived  an  improved  diameter  for  the  asteroid 
Vesta  of  550  ±  23  km,  in  excellent  agreement  with  other 
measurements.  We  have  also  presented  a  measurement 
of  Pallas  of  673  ±  55  km,  slightly  larger  than  previous 
measurements.  However,  Pallas  has  an  angular  extent 
close  to  the  telescope  seeing  limit,  and  the  diameter  may 
be  affected  by  several  systematic  effects  resulting  from 
the  use  of  a  photographic  recording  system.  New  mea¬ 
surements  arc  planned  utilizing  linear  digital  detector 
systems  to  reduce  the  effects  of  these  errors. 

We  have  also  pointed  out  in  this  paper  the  necessity 
of  normalizing  and  centering  each  speckle  frame  before 
reduction,  for  otherwise  the  diameter  derived  may  be 
spuriously  large. 
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